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ABSTRACT 

The Intcrfcrometric Bidimcnsional Spectrometer (IBIS) installed at the Dunn Solar Telescope 
of the NSO/SP is used to investigate the morphology and dynamics of the lower chromosphere and 
the virtually non-magnetic fluctosphere below. The study addresses in particular the structure 
of magnetic elements that extend into these layers. We choose different quiet Sun regions in and 
outside coronal holes. In inter-network regions with no significant magnetic flux contributions 
above the detection limit of IBIS, we find intensity structures with the characteristics of a shock 
wave pattern. The magnetic flux elements in the network are long lived and seem to resemble the 
spatially extended counterparts to the underlying photospheric magnetic elements. We suggest a 
modification to common methods to derive the line-of-sight magnetic field strength and explain 
some of the difficulties in deriving the magnetic field vector from observations of the fluctosphere. 

Subject headings: Sun - spectro-polarimetry, Sun - chromospfiere, Sun - dynamics, Sun - morphology 



1. Introduction 

Understanding the chromosphere and its mag- 
netic structure is important to progress in un- 
derstanding the solar atmosphere as a cou- 
pled phenomenon from the photosphere to the 
corona. A comprehensive physical model of 
the chromosphere must be capable of explain- 
ing a multitude of observations including, e.g.. 
Calcium K grains (Beckers 1968f). the magnetic 
canopy ( Jones and Giovanellil 1983r ) . and chromo- 



spheric oscillation s (jKneer and von Uexkul]|[l995 
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Lites et al. Constructing such a model is 

complicated by many things, and detailed obser- 
vations such as those presented in this work are 
important to constraint this effort. The few avail- 
able chromospheric diagnostics are unfortunately 
subject to complicated formation processes, for 
which - in most cases - the simplifying assump- 
tion of local thermodynamic equilibrium cannot be 
made, causing severe difficulties for the quantita- 
tive comparison between models and observations. 
Furthermore, fluctuations on small spatial and 
temporal scales are very challenging for both mod- 
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Fig. 1. — a) The Ca II infrared line (854.2 nni). Displayed as * are the sampled wavelength points (30) 
of data set 2. the O indicate those of data set 1. The shaded area indicates the wavelength and intensity 
regime within which bisectors for the use with the hybrid bisector-COG were calculated (see text), b) Stokes 
profiles of the strong magnetic features located at [11,16] arcs in Fig. [3]c) (black solid: Stokes V/I; dashed: 
Stokes Q/I; dotted: Stokes U/I). 



eling and observation. Advances in both fields has 
lately led to an increase of attention paid to the 
solar chromosphere, sparking a controversy about 
its "true" physical nature. A rece nt summary of 
this c ontroversy has been given by IVecchio et al. 
(|2009[ ). 

Recent high-resolution observations suggest 
that in the quiet Sun there exists a weak-field 



domain below the classical canopy (jWoger et al 
20061 ). The prominent magnetic fields of the 
canopy are major constituents of the chromo- 
sphere in the stricter sense as it is seen in 
the Ha line. The weak- field domain below 
referred to as "fluctosphere" hereafter 



term first used by IWedemever-Bohm and Woger 
(120081) and previously nain ed "clapotisphere' 



uapotisp 

()Rutten and Uitenbroeklll991l : lRuttCTlll995( l. The 
fluctosphere in numerical radiation hydrodynamic 
models is generated by interference of (acoustic) 
shock waves, which are excited in the photosphere 
and propagate upwards into the layer above (cf. 
Carlsson and Steinlll994[ ). 

At heights around 1000 km above t = 1 
this results in an apparent temperature pattern 
that has similarity to reversed granulation, yet 
changes with typical tim escales between 20-30 s 
( Wedemeyer et al.1l20Q4l ). This pattern is the re- 



sult of the steep temperature increase at the ver- 
tices of the interfering wavefronts, whereas the 
enclosed post-shock regions are cool on average. 

The recent ability to synthesi ze spectra from 

mode ls depicting the fluctosphere (|Wedemever-B6hm et al 
20091 ) allow their direct comparison to obser- 



vations of this atmospheric region, under both 
morphological as well as dynamical aspects. 
Acoustic shocks hav e been analyzed in detail by 
Vecchio et al. I (|2009l) . who use the Ca II infrared 



line at 854.2 nm as diagnostic of what they call 
the "mid-chromospheric" region of the solar at- 
mosph ere. It has been shown by ICauzzi et al 



(|2008[ ) that the Ca II infrared triplet provides a 
convenient and accessible d i agnos ti c of the chro- 
mosph ere. As lJudge et al.l pOOlf ). IVecchio et al 



(|2009[ ) notice a strong influence of magnetic field 
on acoustic processes within that regime of the so- 
lar atmosphere, and conclude that this influence 
may be larger than generally expected. Never- 
theless, they also confirm the existence of regions 
with shock behavior that is similar to what is pre- 
dicted by numerical models with no or at least 
weak magnetic field only. 

A fundamental knowledge of the chromospheric 
magnetic field topology and its dynamics thus 
seems to be the next unavoidable step towards 



2 




^ 4 « 

sat [arvs] 



2 4 G 
9lze lata,] 



i a 5 



2 4 6 




Fig. 2. — Sequence of subframes of a weak-field region in the Ca II infrared line core, in comparison to the 
synthesized model that was degraded with the point spread function of the Dunn Solar Telescope (right 
panel). Intensity is scaled to the same range in all panels. 



a comprehensive understanding of the solar at- 
mosphere. However, the detailed magnetic prop- 
erties of the "mid-chromosphere" are widely un- 
known so far because of the obstacles that need to 
be overcome in observational techniques. This 
in particular applies to the weak-field regions 
in the quiet Sun. It has been suggested that 
magneto-acoustic heating pl ays a role in the 



chromospheric heating (e.g . Hasan et al. 2005t 



Hasan and van Ballegooiien 2008[ ). even though 



the predominant mechanism is still debated. 

Here we present and interpret for the first time 
high spatially resolved spectro-polarimetric mea- 
surements in the fluctosphere above strong pho- 
tospheric magnetic features. We will analyze the 
fine-structure of the line wing and core intensity in 
the Ca II infrared line in comparison to synthetic 
maps from numerical models. The paper is struc- 
tured as follows. In Sect. [2l we present the data 
observed along with the methods used for calibra- 
tion. Section [3] describes the algorithm and result 
of the analysis of the three-dimensional topology 
inherent to a magnetic structure found in our data, 
whereas in Sect. |4]we show the results of our time 
sequence analysis. We conclude with a summary 
of our results in Sect. [5l 

2. Data 

2.1. Data sets 

Data set 1: We observed a quiet Sun region lo- 
cated in a coronal hole at disc-center on May 26, 
2007 using the Call infrared line at 854.2 nm. Sev- 
eral persistant G-band bright points are visible in 
the field of view, forming a network element. The 
overall observation lasted for 26.5 minutes. 



The IBIS narrowband channel was set up in 
spectro-polarimetric mode to scan the line core of 
the Ca II line in two spatial dimensions using six 
modulation states and an exposure time of 35 ms. 
The line was scanned with 17 wavelength steps us- 
ing a step width in wavelength of 4.3 pm with a full 
width at half maximum (FWHM) transmission of 
4.6 pm (see Fig. [T]). This procedure resulted in a 
overall cadence of approximately 27 s for one data 
cube containing two spatial dimensions and a third 
spectral dimension. 

A subregion of the data set displayed for the Ca 
II infrared line core is shown in Fig. [2] for a small 
part of the time sequence. Fig. [3] gives an overview 
over the whole field of view, where total circular 
polarization was defined as 



d\^{v{\)/i{\)Y 



(1) 



and total linear polarization was defined as 



Tc = dAV(Q2(A) + C/2(A))//2(A). (2) 

This data set was observed to acquire further 
insight into the dynamics of the chromosphere 
at locations of strong and weak magnetic flux. 
Unfortunately, after reduction of the spectro- 
polarimatric data of data set 1, the Stokes U com- 
ponent has been found to contain fringes. The 
fringes were compensated, however, at the cost 
of a reduced signal to noise ratio. In addition, 
the Stokes Q component was plagued by a weak 
reflection in the field of view, which hampered 
the Stokes Q signal in the right half of the FOV. 
Overall, the Stokes Q and U components were 
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Fig. 3. — Data set 1: Temporal averages over 26.5 min of a) G-band intensity, b) Ca II line core intensity, c) 
total circular polarisation. The inscribed solid box in a) - c) indicates the location of the subfield presented 
in Fig. [51 the dashed line in therein indicates the position of the simulated slit in Fig. \7\ d) Total linear 
polarisation of the right half of the field of view indicated with the dotted box in a) - c). e) Stokes Q/I 
profiles of enhanced total linear polarization (> 1.9%) in the square inscribed in d). The average of the 
displayed Stokes Q/I profiles leads to the Stokes Q/I profile visible in Fig. [T] 



barely above the noise level, even after spatial 
and temporal averaging. 

Data set 2: Complementary to data within a 
coronal hole, we have observed the Sun on disk- 
center outside of a coronal hole. These observa- 
tions were carried out on September 22, 2008. In 
contrast to the coronal hole observations, no time 
sequence was acquired. However, the exposure 
time was increased to 90 ms for enhanced sensi- 
tivity. In addition to the Ca II infrared spectra, 
the Fe I line at 630.2 nm was observed with IBIS 
subsequently with a time difference of 1 minute. 
For these observations, the number of wavelength 
points was increased to 30 in the Ca II infrared line 
to cover a larger part of the spectrum (Fig.[T|). The 
Fe I line was scanned with 30 wavelength points 



with a stepwidth of 3.2 pm. The observations thus 
includes both the 630.1 and 630.2 lines, as well as 
the telluric blend spearating these two lines. The 
FWHM transmission of IBIS at this wavelength is 
about 2.3 pm. 

This data set was acquired complementary to 
data set 1 to gather additional information, which 
is used to check the magnetic field information re- 
trieved from the Ca II infrared line. 

2.2. Data calibration 

The IBIS instrument consists of two channels, a 
narrowband channel with two tunable Fabry-Perot 
interferometers and a broadband channel for ref- 
erence (jCavallini 20061) . The broadband images 
of IBIS were calibrated with the standard proce- 
dure of dark subtraction and gain table application 
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and reconstructed using the speckle image recon- 
struction algorithm modified for hig h-order adap- 
tive op tics corrected data detailed in IWoger et alJ 
()2008l) . 



The IBIS narrowband data were calibrated us- 
ing the following procedure. In a first step, the 
gathered flat and dark images were averaged. Be- 
cause the Fabry-Perot interfcrometers of IBIS are 
located in a coUimated beam, each pixel in a sin- 
gle narrowband exposure is affected by a blueshift 
( Cavallini 20061 ). The magnitude of the blueshift 
was estimated from the dark corrected data cube 
of average flat images by determining for each spa- 
tial pixel the line core position to sub-pixel accu- 
racy. The resultant two-dimensional blueshift map 
then fitted with its theoretical form, a paraboloid. 
In a subsequent step, the mean line profile was 
computed by averaging the blueshift corrected sin- 
gle profiles. To compute the gain table, each line 
profile at a pixel position in the data cube was 
divided by the mean profile which was shifted to 
match the line core position at that specific pixel. 
After the application of the gain table to the data, 
the narrowband images of IBIS were destretched 
using the broadband images in combination with 
their reconstruction. Finally, all line profiles were 
blueshift corrected, leading to a pixel artefact free 
data cube, which served as input for the Stokes 
parameter calibration. 

The polarization calibration is accomplished us- 
ing procedures originally developed f or the Ad 



vanced Stokes Polarimeter (ASP) (e.g. iLites et al 



1993af ). that were adapted for usage with IBIS 
spectro-polarimetric data. The procedure is wave- 
length independe nt, and has be e n pub lished for 



IBIS in detail by IViticchie et al 
larization calibration includes 



I 2009[) . The po- 
removal of resid- 



ual cr osstalk between the Sto kes parameters. Re- 
cently, [judgT^ni^eLS ( 2009[ ) has compared IBIS 
spectro-polarimetric data simulatenously observed 
to that of the spectrograph (SOT/SP) onboard 
the HINODE satellite, the result of which indi- 
cates that the calibration processes yield very sim- 
ilar results. A publication describing in detail the 
complete calibration of spectroscopic and spetro- 
polarimctric da ta acquired with IBIS is c urrently 



3. Topology of magnetic network elements 

For data set 2, we derive the topology of strong 
magnetic fields reaching from the photosphere into 
the chromosphere. An example of the data quality 
is given in Fig. [4l A strong magnetic feature con- 
sisting of photospheric bright points can be seen 
in the intensity map at a continuum wavelength 
of 828.04 nm. The chromospheric counterpart is 
highly structured and consists of many fibrils. Us- 
i ng th e method put forward by Rees and Seme! 
(|l979l ). and later also lUitenbroekl pOOSh . we derive 
the LOS velocity and the magnetic flux density. 
The results for the 3D topology of the magnetic el- 
ement are very similar for the Fe I and the Ca II IR 
line (see Fig. |4l panels d) and h)). It appears that 
both predominantly map the magnetic flux in the 
photosphere. This is due to the fact that the result 
of the COG method represents an average over the 
entire Ca II IR line profile which encodes the chro- 
mospere only within a very small wavelength in- 
terval around the line core, whereas the line wings 
- that encode mostly the photosphere - make out 
the major part of the profile. Thus, the result is 
significantly biased towards lower atmospheric re- 
gions. 

3.1. The hybrid bisector-COG method 



The center of gravity (COG) method (jUitenbroek 



in preparation (|Tritschler and et al. Illooi). 



20031 ) to derive LOS magnetic flux is very robust, 
but has the disadvantage that height fluctuations 
in the field cannot be recovered. For such an anal- 
ysis, usually a bisector analysis is employed. The 
latter is unfortunately sensitive to noise. In order 
to increase the reliablility of our topology analysis, 
we devise a hybrid method combining the bisector 
method with the COG method. 

In the algorithm, the following steps are per- 
formed: 

1. Compute the wavelengths at which the red 
and blue line wing of both the I-l-V and I-V 
spectra have the same intensity. This can be 
pictured as finding the two points (Aa and 
Xb) of intersection of the line profile with a 
horizontal line. 

2. Use the two acquired wavelength ranges to 
compute the COG of the I-f-V and I-V spec- 
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Fig. 4. — Data set 2, acquired on September 22, 2008. a) IBIS continuum image, b) Fe I (630.2 nm) core 
image, c) line-of-sight velocity of Fe I d) linc-of-sight magnetic flux density in Fe I, e) Ca II (854.2 nm) core 
intensity, with a contrast enhanced close-up of the marked fibrils in the subpaneL f ) Ca II wing intensity (+ 
0.073 nm), g) line-of-sight velocity of Ca II and h) line-of-sight magnetic flux density in Ca II. Ve locities and 
magne tic flux densities have been calculated using the genuine COG method analyzed in detail in lUitenbroekl 
(l2003h . 
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tra with 



A+ = 



/;;dAA{(I+V)(A,)-(I+V)(A)} 
/,';dA{(I+V)(A,)-(I+V)(A)} 



(3) 



for I+V and analogously A_ for I-V. 



3. The value for the magnetic flux density along 
the line-of-sight is then computed from the 
difference of the COG positions of the I+V 
and I-V profiles, using 



A, 



A_ A TT m c 



\ 2 ' 

e gL Xq 



(4) 



where Aq is the central wavelength of the 
line, A± are the wavelength positions of 
the centroids of I±V, gi is the effective 
Lande factor, and e and m are the elec- 
tron char ge and mass, res pectively, in SI 
units (see lUitenbroekl |2003| . and references 
therein) . 

Step 1 is needed for two reasons: First, the selec- 
tion of many different intensity values allows to 
analyze the magnetic flux density at different for- 
mation height ranges, which is also true for the 
genuine bisector method. Secondly, an uneven cut 
of the I+V and I-V would bias the measurement of 
the COG value. However, it is obvious that this 
method only allows an analysis of the magnetic 
flux density in non-emitting lines, as an emission 
peak would make Step 1 and Step 2 very diffi- 
cult. An emission peak would create four instead 
of two points of intersection and also bias the com- 
putation of the COG. Another weakness of the 
described method is sensitivity to asymmetries in 
the Stokes V profiles. Such asymmetries are found 
e.g. in a detaile d analysis of a Quiet Sun region 
near network by IPietarila et al. ( 2007 ) on a com- 
mon basis, and can lead to biased measurements 
of magnetic flux. In what follows, we test the pro- 
posed method with synthetic profiles that inhibit 
similar behavior. 

3.2. Check of method 

We have verified the feasibility of this approach 
by applying it to a magnetohydrostatic model of 
a funnel-like expanding fluxtube. It originates in 
the deep photosphere with a footpoint of 500 km 
in diameter and has a magnetic field strength of 



2500 Gauss. NLTE spectra for all four Stokes pa- 
rameters were subsequently synthesized for the 
Ca II IR line. The result is shown in Fig. [H The 
absolute magnetic flux density was accurately re- 
covered where magnetic field was present, as was 
the photospheric topology of the fluxtube. No sig- 
nificant jumps occur between the results computed 
from adjacent intensity levels. This is likely due 
to the COG algorithm that reduces noise in the 
measurement by including many more points in 
the derivation of the magnetic flux density than 
the simple bisector method. Interpretation of the 
three-dimensional topology in the higher layers be- 
comes difiicult due to the fact that a bisector does 
not necessarily correspond to a thin atmospheric 
layer. 

The effective formation height was derived for 
the Ca II infrared line by averaging the geometri- 
cal heights of the locations where t\ ^ 1 over each 
of the wavelength intervals computed in Step 1. 
In general, the formation height increases from 
line wing towards the line core. In addition, the 
"Wilson effect" causes the height of ta = 1 to 
differ in magnetic and non-magnetic regions. The 
strong magnetic fiux causes the opacity within the 
fiux tube to be lower compared to the ambient 
medium, resulting in contributions to line core sig- 
nal from layers at about 200 km lower compared to 
the T\ = I height in the non-magnetic surround- 
ing. The Wilson effect has been compensated in 
Fig. [5] b) , which demonstrates the well recovered 
topology of the flux tube. 

The comparison of the input model in panel a) 
and the reconstruction in panel b) of [5] serves now 
as reliability test and error estimate for the pro- 
posed method. 

The recoverable height interval is here re- 
stricted by the line properties such as forma- 
tion height and the chosen wavelength intervals 
in Step 1. Qualitatively, the topology of the mod- 
eled fluxtube was recovered. The regions without 
magnetic field in the input model show a weak 
signal in the reconstruction corresponding to field 
strengths of up to (6 ± 1) G. These artefacts can 
primarily be attributed to the fact that the Stokes 
V spectra exhibit features like spikes when dis- 
continuities occur along the LOS, biasing the I+V 
and I-V spectra and introducing a fake magnetic 
signal. 
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Fig. 5. — Result of the accuracy test of the proposed hybrid bisector-COG method, a) input model, and 
b) recovered flux density. The height scale in b) was derived using calculations of ta = 1. The gray scale 
encodes the magnetic flux density and is the same for both panels. 



3.3. Reconstructed field topology in data 
set 2 

We now analyze dataset 2 and use 101 intensity 
values ranging from 18.5% (Ca II IR line core) to 
32% of the local continuum intensity (Fig. [1]). We 
used a median filter with 5x5 pixels (correspond- 
ing to 0.85x0.85arcs^) to reduce salt and pepper 
noise introduced by the measured spectra. 

First we apply the COG method to the Fel 
line data and obtain a magnetic flux density of 
about 350 Gauss for the strong feature in the FOV, 
likely because the small-scale magnetic elements 
have not been resolved, yielding a decreased field 
strength because the signal is averaged over the 
area of a pixel. The COG method is sufficient here 
because this spectral line is formed over a compar- 
atively small height range in the atmosphere. The 
formation of the Call IR line, on the other hand, 
extends over more than a 1000 km. 

Nevertheless, the computed magnetic flux den- 
sity can serve as a well-tested value for those re- 
trieved using the Ca II IR line. Indeed, the re- 
sult fits well with the 250 Gauss obtained from the 
measured Ca II IR line using the same method, 
which started at ~50% of continuum intensity; the 
formation height of the magnetic signal of the Ca 
II IR line cxpcctedly appears to be higher than 
that of Fc I. The value drops further to 125 Gauss 
with our hybrid-COG method when using bisec- 
tors starting at 32% of the local continuum in- 
tensity. This is not surprising because we have 
restricted ourselves even further to the line core. 



The result of the hybrid bisector-COG analysis of 
our data can be seen in Fig. [51 

It is worth to point out two situations visible 
in the data. A patch of magnetic flux directed in 
the direction of the observer seems to disappear 
when restricting oneself to the line core with the 
derivation (Fig. [HI location A), while at the same 
time, a filament like structure appears (Fig. [51 lo- 
cation B) . Obviously, the patch of location A does 
not reach far up into the solar atmosphere, likely 
connecting to regions with opposite polarity in the 
photosphere. The filament like structure, on the 
other hand, appears in a location of enhanced core 
intensity of the Ca II IR line (see also Fig. [4]). To 
our knowledge, this is the first time a filament-like 
structure has been detected in magnetic signal. 

The confined photospheric magnetic structure 
appears to become fragmented at higher atmo- 
spheric levels. Overall, it appears that the struc- 
ture of the magnetic flux density along the LOS, 
i.e., essentially the vertical component, becomes 
smaller and weaker. This can be explained with 
a magnetic funnel expanding with height. In that 
case, the magnetic field becomes increasingly hor- 
izontal with height, so that the vertical flux com- 
ponent decreases while the horizontal component 
increases. 

Due to instrumental limitations, weak horizon- 
tal fields cannot be detected in data set 2. The 
failure to measure the horizontal field in total lin- 
ear polarization can be explained with the finite 
sensitivity of these particular IBIS observations - 
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Fig. 6. — Cuts through the recovered 3D topology of the magnetic flux density by applying the hybrid 
bisector-COG method to the Ca II IR line scan of dataset 2. Panel a) is at 32% of the local contimmm 
intensity (photospheric) (see Fig. [T]), b) at 20.2% and c) at 18.5% (fluctospheric), with a contrast enhanced 
closc-up of the fibris marked by the box in a subpanel. At location A, a patch of magnetic flux along the 
line-of-sight seems to disappear with increasing height (also see subpanel). Location B is an example for the 
appearance of a filament like structure with height. 



no temporal averaging was possible as data set 2 
only consists of a single line scan. 

Using temporal averaging it is possible to in- 
crease the signal to noise ratio in the spectro- 
polarimetric data by decreasing photon noise, 
while losing the ability to perform any dynamic 
analysis. To detect horizontal magnetic fields, 
we have temporally averaged data set 1. The 
integrated exposure time of the resulting maps 
is equivalent to 12.5 s. The result is shown in 
Fig.[3]a)~d). The Stokes parameters averaged over 
the time sequence are displayed in Fig.[T]b). There 
are strong magnetic features (around x= 11 arcs, 
y= 15 arcs), which show up as distinct signals in 
the total circular polarization signal with a di- 
ameter of 1 arcs. This feature can also be seen 
at the same location in the photosphere, partic- 
ularly well in G-band images. Only the summed 



up data have sufficient signal to noise ratio to 
study the change in topology of the magnetic flux 
density that extends from the photosphere into 
the chromosphere. The photospheric footpoints 
of the magnetic features are advected with the 
convective motions. As the duration of the time 
sequence is rather short, the features roughly re- 
main at the same location. Thus, in the averaged 
G-band images the location of the magnetic fea- 
tures are clearly visible (Fig. [3] a)). Temporal 
averaging reveals a structure in total linear polar- 
ization encirceling the magnetic elements visible in 
total circular polarization (Fig.[3]d)). This struc- 
ture supports the view of a horizontal magnetic 
flux component that increases from photosphere 
to chromosphere: as the magnetic flux becomes 
weaker with increasing height and distance from 
its footpoint in the photosphere, it becomes un- 
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measurable even with an equivalent exposure time 
of 12.5s. 

3.4. Discussion 

When analyzing the expected polarimetric sig- 
nal of e.g. the hydrostatic model of an expand- 
ing flux tube described above and displayed in 
Fig. [5] a), a very weak total linear polarization 
signal is expected (of the order of 0.03% of the 
measured intensity) , which is below the sensitivity 
of IBIS even with 12.5 s integrated exposure time. 
The measured linear polarization signal is of the 
order of 0.1% of the measured intensity (Fig. [T]), 
which is barely above the noise level, whereas am- 
plitude and spatial extent of the input model flux- 
tube at photospheric heights flt well with the mea- 
sured total circular polarization signal in our data 
sets. 

Possible explanations for this finding are sum- 
marized in what follows. 

• The origin of the measured signal cannot be 
attributed easily to a certain height layer in 
the solar atmosphere. Significant contribu- 
tions to the polarization signal may originate 
from atmospheric regions with (a) strong 
magnetic field for which the contribution 
function is small, and (b) weak magnetic 
field for which the function has large val- 
ues. The total circular polarization signal 
may therefore originate from a different at- 
mospheric layer than the total linear polar- 
ization signal. These issues may be the rea- 
son for the agreement of the total circular 
polarization signal with the model of an ex- 
panding flux tube, as its signature originates 
from deeper layers than the fluctosphere. In 
these regions the model may predict the flux 
tube expansion correctly. 

• The "large" values of the total linear po- 
larization may indicate that loop-like struc- 
tures prevail which connect to magnetic el- 
ements in the photosphere. In this case 
the total linear polarization signal (and thus 
the horizontal magnetic field) may originate 
from atmospheric regions that lie well be- 
low the flu ctosphere in the u p per photo- 
sphere (e.g 'Steiner et al.' '2008'; ' Lites et al 
[2OO8 : Ishikawa and Tsuneta 2009;) In those 
regions, the magnetic field lines are more 



confined and produce a stronger polarization 
signal. 

The model of a static fiux tube that expands 
as a funnel-like structure into the chromo- 
sphere likely is too simplified to create an ac- 
curate picture of the dynamic chromosphere. 
Recent observations and simulations sug- 
gest that magnetic field elements may have 
a more complex and entangled structure 
than thought before (see recent reviews by 
Wedemever-Bohm et al.ll2009t Vecchio et al 



2009 : Pietarila et al.l 20091 and references 



therein). 

4. Small-scale intensity pattern in inter- 
network regions 

The short intensity map sequence displayed in 
Fig. [2] is characterized by very short-lived grain- 
like features of enhanced intensity. Their sizes 
are of the order of an arcsec or less. While such 
Ca grains are well-known phenomenon, the sub- 
structure of the background has been less studied 
- due to instrumental limitations. It is the com- 
bination of a very narrow transmitted wavelength 
range, a short exposure time, and the high spa- 
tial resolution due to the well performing adap- 
tive optics system that enable us resolve the pat- 
tern between the grains. We now clearly see dim 
strands connecting the grains. While the grains 
appeared as isolated features before, they now ap- 
pear to be the nodes of a mesh-like pattern with 
rather dark regions in-between. Short-lived "cal- 
cium grains" have been known for a long time (e.g. 
Kneer and von Uexkulll 119931 ) not only in Ca II 
K, but also in the Ca II IR line. The formation 
of "calcium gr ains" has been explained using ID 
simulations bv lCarlsson and Stein (1997) as prop- 
agating shock waves. Indeed we find grains in con- 
nection with a clear signature of shock waves in the 
time variation of the spectrum (Fig. [7|) . This sig- 
nature is typical throughout the weak-field regions 
of data set 1. 

The intensity pattern at small spatial scales 
resembles the pattern produced by interaction 
of propagating shock waves as seen in radiation 



hydro dynamic simulations of IWedemever et al 
(l2004[ ). We therefore compare the observed line 



core maps with the synthetic inten s ity m aps de- 
scribed in IWedemever-Bohm et al] (|2009[) . The 
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latter were calculated wit h the non-LTE radia- 
tive transfer code MULTI ( Carlsson 1986h from a 
time-dependent radiat ion hydrodynamics s imula- 



tion with CO^BOLD (|Freytag et all 120081 ). The 



tern, which is also connected to a gas temperature 
distribution with a hot shock wave component and 
a cool background. Already the simulations by 
Skartlien et all (|2000f ) produced a shock-induced 



line core intensity originates from a height range of 
around 1000-1500 km (see the contribution func- 
tions by lUitenbroeklbood ). The subregion shown 
in Fig. [2] was chosen to match the size of the 
model. The synthetic map in the figure was de- 
graded by convolution with a theoretical point 
spread function (PSF) of the Dunn Solar Tele- 
scope . As described in IWedemever-Bohm et al.l 
(|2009l ). this PSF is rather optimistic. It is hkely 
that the instrumental straylight is underestimated 
and that residual atmospheric distortions are still 
present even though the observations were ac- 
quired with the high-order adaptive optics system 
operating near the diffraction limit, leading to a 
too high contrast compared to the observations. 
And still most of the filigree fine-structure in the 
original synthetic maps is lost after application 
of the PSF. What remains is very similar to the 
observation: A mesh-like pattern enclosing dark 
regions and bright grain-like features at the ver- 
tices of the mesh. In the model, the brightening 
occurs at locations with high gas temperature as a 
result of the compression of the gas in the collision 
zone between neighboring shock waves. 

The aforementioned simulations do not take 
into account magnetic fields and thus represent 
highly idealized conditions. The observations dis- 
played in Fig. [2] were taken in a very quiet re- 
gion within a coronal hole and are therefore most 
likely as quiet as they could possibly be. In- 
deed, the line core images exhibit no fine-structure 
with "fibrils" or "threads" which are otherwise 
seen in the vicinity o f stronger magnetic fields (cf. 
data se t 2 and, e.g.. iRouppe van der Voort et al 



20051 : IVecchio et al.ll2009l ). Although no clear pat- 
tern is discernible neither in the total circular 
nor in the total linear polarization measured in 
the subregion (see Fig. [3]), weak magnetic fields 
at or below the detection limit could well be an 
integral component of the fluctosphere even un- 
der su£h_extreme_conditio^ The simulations 
by ISchaffenberger et~ ( 2006[ ). which include a 
weak magnetic field in the fluctosphere, still ex- 
hibit a shock wave pattern similar to those seen 
in the field- free simulations. The simulations by 
Martinez-Svkora et al. ( 2008[ ) show a similar pat- 



pattern can, although the relatively coarse grid 
spacing makes it harder to discern. 

In Fig. [7] we compare the spectra of the non- 
magnetic model degraded with the point spread 
function of the Dunn Solar Telescope with one of 
our measurements. The figure shows the tempo- 
ral evolution of the spectra in the magnetic ele- 
ment visible at [11,16] arcs, e.g., in Fig. [3]c) and 
at [24,3] arcs, a region with no appearant magnetic 
signal (Fig. [7]b). The spectrum from the region 
without appearant m agnetic flux is almost identi- 
cal to those visible in Vecchio et al. ( 20091 ) , show- 
ing frequent and regular shock events. The spec- 
trum in the magnetic element (Fig. [7] a) shows a 
different picture: here, shocks happen much rarer 
and irregular. It demonstrates the influence of 
magnetic fields on the generation and propagation 
of s hock waves as it has been frequently observed 
(cf. I Vecchio et al]l2009^ ■ In addition, events in 
the spectrum appear to be mostly simultaneous in 
blue and red wing, whereas changes in the non- 
magnetic temporal spectra (Fig.[7]b) happen with 
temporal delay indicating acoustic events. 

Shock events (e.g., at < = 530s) produce sim- 
ilar structures in synthetic and measured spec- 
tra. Both, the spatial and spectral extension of 
the short lived brightening at a position of a; = 
2.5 — 3 arcs in the measured and at x = 6 arcs 
in the synthetic spectrum are very similar. In the 
synthetic spectrum, a brightening of about 0.5 arcs 
in size appears in the line core. 

Owing to the very steep change in effective for- 
mation height range with wavelength from core to 
wing, the fluctospheric pattern is only visible in 
the line core maps. Already a few wavelength posi- 
tions further into the wing, the intensity pattern is 
dominated by contribution from the layers below. 
The imprint of reversed granulation thus already 
appears when observing at wavelengths close to 
the line center (AA « ±30 pm. The rather high 
velocities in the fluctosphere cause Doppler shifts 
that in principle mix the intensity signature from 
different layers from pixel to pixel and possibly on 
a sub-pixel scale. One should thus beware that 
despite the narrow transmission range of the IBIS 
instnmient the Ca II 854.2 nm line core maps could 
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Fig. 7. — a) Time variation of the spectrum of the Ca II IR hne at the location [11,16] arcs in the field of view 
of dataset 1 (magnetic element), b) Like panel a), except for a weak-field location at [24,3] arcs, c) Spatial 
spectrum of dataset 1 at t=530s along the slit indicated in Fig. [3]b) with a length of ~7.7arcs. d) Spatial 
variation of the synthesized model degraded using the point spread function of the Dunn Solar Telescope 
along a slit at position a; = 2.75 arcs, with a slit length of 7.7 arcs. 



contain significant contributions from an extended 
height range in the atmosphere. Nevertheless, we 
are confident that the intensity pattern observed 
here is indeed a fluctospheric shock wave pattern 
as seen in numerical simulations and not the re- 
versed granulation pattern at mid-photospheric 
heights. Already visual inspection of the image 
sequence indicates that the two-dimensional pat- 
tern evolves faster than known for reversed gran- 
ulation. 

This finding is quantified here in terms of the 
typical time scales with which the two-dimensional 
pattern seen at different wavelength positions in 
the Call infrared line evolve. We calculate the 
two-dimensional auto-correlation function 

Af{t) : ^ 



■lit 

Nt~t-1 

E 

if=0 



T,uLo^{{F{ix,iy,it)f)x,y 

{{F{i,,iy,it)-F) X (5) 
(Pi 

it + r) - F)) 
for the discrete, time-dependent function F(x, y, t) 



sampled with Nt time steps, which here is an in- 
tensity map sequence at a fixed wavelength. The 
brackets ( • )x,y represent the average over the spa- 
tial elements and F ~ {F{x,y,t))x.y.t represents 
the average of all elements. It should be noted, 
that two images with time lag <ovoi are identical 
if ^F(icvoi) = 1 or perfectly anti-correlated for 
■AF{tcvo\) = — 1- In case -4i?(iovoi) = 0, two im- 
ages with time lag ievoi show no correlation at all 
and the scene has changed completely so that no 
structure related to the first image is present any- 
more. The decay time scale tcvoi of the pattern 
is commonly defined as the point at which Af{t) 
has dropped to AF{tcvo\) = 1/e- The values of 
ievoi were computed for continuum, line wing and 
line core intensity time sequences of the Ca II in- 
frared line in weak-field regions. 

For the Call IR line core at 854.216 nm, the 
correlation time is here evaluated to ip^™ """^ ~ 
59 s. It matches the number found in previ- 
ous observations of the fluctosphere using the 
Call K li ne at 393.3 nm with a broader FWHM 
of 30 pm (jWoger et all l2006l ). In that work, a 
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rapidly changing pattern governs the correlation 
time which is very similar to the pattern detected 
in this work in the Call IR line. In contrast, 
the intensity pattern observed in the line wing 
of the Ca II IR line at 854.246 nm in data set 1 
clearly exhibits reversed granulation. That pat- 
tern evolves on a scale of » 2.5 min. 

iLeenaarts and Wedemever-Bohml ( 20051 ) used the 



observed intensity in the Ca II H line (396.883 nm) 
to analyze the correlation times of the 2D struc- 
tures visible. Their observations - obtained using 
an interference filter with a full width at half max- 
imum (FWHM) of 0.1758nm - show mainly the 
structures of reversed granulation with t'^^^i ~ 
2 minutes. This is caused by the rather broad 
filter range, which includes a significant part of 
the line wings. Consequently, the contributions 
of the brighter wings to the integrated intensity 
is much larger than for the comparatively dark 
core. This way the reversed granulation pattern, 
which is formed in the middle photosphere, is 
the dominant pattern in the observed intensity 
maps. Nevertheless, because the filter was cen- 
tered on the line core of the Ca II H line, contri- 
buti ons from the fast changing fluctospheric lay- 



ers (jWoger et al 



2006 ) are mixed with those from 



the lower layers ( Uitenbroek 20061 ). This leads to 
shorter correlation times than those of pure re- 
versed granulation. Separating the faster evolving 
fluctospheric pattern from reversed granulation re- 
quires filters with a very narrow transmitted wave- 
length range, like the FWHM of 4.6 pm for IBIS 
as used here. 

Figure [8] a) shows a scatter plot of evolution 
time scale versus total circular polarization. The 
correlation function was computed in fields with 
3.74^ arcs^ field size to sample our FOV on a finer 
grid. On average, the resultant evolution time 
scales were not changed significantly by the de- 
crease of the subfield size. We now discriminate 
between different types of subfields by imposing 
thresholds for evolution time scale and total circu- 
lar polarization. The individual subfields are cor- 
respondingly classified. The different gray scales 
in Fig. [8] b) show the distribution of the differ- 
ent types in the FOV. The subfields with the 
longest correlation time (white patches) arc con- 
nected to small-scale, long lasting structures that 
are hkely related to n etwork magnetic field (see 
also lWoger et al.ll2006l ). 



Not only vertical fields play a role in the ap- 
pearance and dynamics of structures visible in the 
Ca II IR line, but also horizontal fields, as the gray 
subfields are primarily located above the struc- 
ture in total linear polarization visible in Fig. [3] 
It has long been known that the magnetic flux 
plays an importan t role in emer ging intensity in 
the chromosphere (iBeckerd 1968 ) , and oscillations 
(jLites et al.l Il993br ). as it impacts e.g. the gas 



pressure and thus the contribution function. The 
structures with short life times prevail in regions 
with magnetic flux below the detection limit of 
IBIS. Obviously, the magnetic field has an impact 
on the dynamics of the measurable intensity struc- 
tures that exist in the "mid-chromosphere" . 

However, the cadence of 27 s may not be suffi- 
cient to reliably detect time scales that are lower 
than twice the cadence of the data set. A shorter 
cadence is certainly desirable, in particular be- 
cause numerical simulations indicate time scales 
that are much shorter. Such observations will be- 
come possible in the near future - an IBIS upgrade 
with faster cameras is currently in progress. 

5. Conclusion 

Spcctro-polarimetric observations of the solar 
atmosphere with a focus on heights equivalent to 
the low chromosphere have been presented in this 
work. We find that fluctospheric regions of the 
solar atmosphere, i.e. the weak field domain be- 
tween the photosphere and the magnetically dom- 
inated chromosphere, display an intensity pattern 
that shows structures of similar spatial scales as re- 
versed granulation but is evolving on much faster 
time scales. From comparison to the predicted 
synthesized intensity of three-dimensional radia- 
tion hydrodynamic models of the fluctosphcre, we 
conclude that the pattern is the intensity signa- 
ture of propagating, interacting shock waves. The 
pattern's dynamic time scale of 59 s is comparable 
to the 20-30 s of temperature cuts at the height 
of 1000 km above Tcont = 1, when taking into ac- 
count that seeing and finite spectral passband of 
the instrument will lead to prolonged correlation 
times. 

We propose a method to derive the 3D topol- 
ogy of the absolute LOS magnetic field fiux per 
pixel. The method works for simple cases in non- 
emitting lines like the Ca II IR line in the Quiet 
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Fig. 8. — a) Scatter plot of correlation time (computed in subfields of 3.74^ arcs^ size) versus mean total 
circular polarization using data set 1. b) Origin of the data points displayed as symbols in a), overlayed 
over the temporal average of the total circular polarization of data set 1: O originate from the dark areas, 
+ originate from the gray subfields and * originate from the white subfields (marked with a white box). 



Sun. For such cases, we were capable to recon- 
struct the topology of a strong magnetic structure 
that extends from photosphere to chromosphere. 
The structure becomes weaker with height, but 
does develop a filament like structure when recov- 
ering the field using the line core. The method is 
not capable to infer filling factors, and thus can- 
not reproduce the field strength of an unresolved 
magnetic element. 

A complementary dataset with higher signal-to- 
noise ratio shows structures in total linear polar- 
ization encirceling strong magnetic elements. This 
strongly suggests a flux tube funnel structure in 
the chromosphere. However, the location of the 
horizontal field remains unclear, and a static flux 
tube funnel appears to be a too simplified model 
for magnetic fields in the chromosphere. 

Highly spatially resolved measurements of the 
magnetic fields in the chromosphere are important 
when trying to gather a deeper understanding of 
the chromospheric energy balance. The work pre- 
sented here clearly shows that, despite the usage 
of advanced techniques, the signal-to-noise ratio 
of currently available data prohibits the detection 
of weak field strengths and definitely not at the 
temporal resolution which is implied by the inten- 
sity variations seen in chromospheric /fluctospheric 
diagnostics such as the Ca IR line core. Ulti- 



mately, telescopes with significantly larger aper- 
tures are needed to get a more detailed picture of 
the chromosphere. The difficulties of the observa- 
tion of weak magnetic fields in the chromosphere 
are only surpassed by the interpretation of such 
data. In general, the height origins of the Stokes 
signals are unclear and do not only depend on con- 
tribution function but also on height distribution 
of the magnetic field strength. In addition, dis- 
continuities and steep gradients in magnetic field 
strength are difficult to detect and remain a funda- 
mental problem of all methods to recover the mag- 
netic field vector from measurements with meth- 
ods available today. 
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